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ABSTRACT 

Supernova (SN) 2006gy was a hydrogen-rich core-coUapse SN that remains one of the most himinous 
optical supernovae ever observed. The total energy budget (> 2x10^^ erg radiated in the optical 
alone) poses many challenges for standard SN theory. We present new ground-based near-infrared 
(NIR) observations of SN 2006gy, as well as a single epoch of Hubble Space Telescope (HST) imaging 
obtained more than two years after the explosion. Our NIR data taken around peak optical emission 
show an evolution that is largely consistent with a cooling blackbody, with tentative evidence for a 
growing NIR excess starting at day ~100. Our late-time Keck adaptive optics (AO) NIR image, taken 
on day 736, shows little change from previous NIR observations taken around day 400. Furthermore, 
the optical TJST observations show a reduced decline rate after day 400, and the SN is bluer on day 825 
than it was at peak. This late-time decline is inconsistent with ^^Co decay, and thus is problematic for 
the various pair-instability SN models used to explain the nature of SN 2006gy. The slow decline of 
the NIR emission can be explained with a light echo, and we confirm that the late-time NIR excess is 
the result of a massive (>10 M©) dusty shell heated by the SN peak luminosity. The late-time optical 
observations require the existence of a scattered light echo, which may be generated by the same dust 
that contributes to the NIR echo. Both the NIR and optical echoes originate in the proximity of 
the progenitor, ^lO'^^ cm for the NIR echo and <20 pc for the optical echo, which provides further 
evidence that the progenitor of SN 2006gy was a very massive star. 
Subject headings: supernovae: general — supernovae: individual(SN 2006gy) 



1. INTRODUCTION 

At the time of discovery, supernova (SN) 2006gy was 
the most luminous SN ever found (Ofek et al. 2007; Smith 
et al. 2007). SN 2006gy generated a great deal of inter- 
est; in addition to being ~100 times more luminous than 
a typical Type II (hydrogen-rich, core-collapse) SN at 
peak, it exhibited a long rise time (~70 day) and slow 
decline, leading to speculation that it may have been the 
first observed example of a pair-instability SN (Ofek et al. 
2007; Smith et al. 2007) or a pulsational pair-instability 
SN (Woosley, Bhnnikov, & Heger 2007). 

SN 2006gy was classified as a Type Iln SN (see Schlegel 
1990 for a definition of the Type Iln subclass and Filip- 
penko 1997 for a review of its spectral properties) based 
on the relatively narrow emission features present in the 
early-time SN spectrum (Ofek et al. 2007; Smith et al. 
2007). Some Type Iln supernovae (SNe Iln) are known 
to be overluminous relative to their typical SN II coun- 
terparts, probably due to the collision of fast-moving 
SN ejecta and a dense, possibly clumpy, circumstcUar 
medium (CSM; Chugai & Danziger 1994). In a com- 
panion paper (Smith et al. 2009b), a detailed spectro- 
scopic comparison of SN 2006gy is made to other SNe Iln. 
SN 2006gy is unique within the SN Iln subclass, how- 
ever, because typical interaction models cannot explain 
its early-time behavior, suggesting the need for alterna- 
tive models for this particular object (Smith et al. 2007; 
Woosley, Bhnnikov, & Heger 2007; Nomoto et al. 2007; 
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Smith et al. 2009b). 

Pair-instability SNe (Barkat, Rakavy, & Sack 1967; 
Rakavy & Shaviv 1967; Bond, Arnett, & Carr 1984) 
are expected to occur in very massive, low-metallicity 
stars, such as those that may have been present in the 
metal- free environment of the very early universe (e.g., 
Abel, Bryan, & Norman 2000). The detection of a pair- 
instability SN in the comparatively local universe, then, 
could potentially reveal a great deal about the first gen- 
eration of stars. The light curves of pair-instability SNe 
are expected to exhibit a relatively slow rise, followed by 
a broad turnover after the peak, and a peak luminos- 
ity that is considerably larger than that of typical SNe 
(Scannapieco et al. 2005). Quahtatively each of these 
characteristics matches those observed for SN 2006gy. 

Other more luminous SNe have been announced since 
the discovery of SN 2006gy: SNe 2005ap (Quimby et 
al. 2007) and 2008es (Miller et al. 2009; Gezari et al. 
2009). This suggests that while these events are rare, 
there may be a new subclass of very luminous supernovae 
(VLSNe). The peak luminosity and photometric evolu- 
tion of both SNe 2005ap and 2008es are difficult to ex- 
plain via the pair-instability model (Quimby et al. 2007; 
Miller et al. 2009), implying that peak luminosities > 
few xlO''^ erg s^^ are possible without a pair- instability 
explosion. There is a wide diversity of alternative mod- 
els that have been developed to explain the early-time 
observations of SN 2006gy. Smith & McCray (2007) ar- 
gue that the peak luminosity and light-curve evolution 
can be explained via the thermalization of shock energy 
deposited into a massive (~10 Mq), optically thick shell. 
Similar explanations were offered for SNe 2005ap and 
2008es (Quimby et al. 2007; Miller et al. 2009). Based 
on a model of the light curve near peak, Agnoletto et al. 
(2009) suggest that the combination of ejecta colliding 
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with dense clumps in the CSM and ^1 3 Mq of '"'^Ni are 
responsible for the early-time luminosity of SN 2006gy. 
Nomoto et al. (2007) are unable to match the light curve 
with a standard pair-instability SN model; however, they 
can reasonably reproduce the (<400 day) light curve of 
SN 2006gy using a pair-instability model where they arti- 
ficially reduce the cjccta mass, such that the radioactive 
heating is less than 100% efficient. Woosley, Blinnikov, 
& Heger (2007) use a model of a pulsational pair insta- 
bility within a massive star to explain the light curve of 
SN 2006gy. 

Some of the models make predictions for the late-time 
behavior of SN 2006gy. Those with a large yield of '"'^Ni 
(Nomoto et al. 2007; Smith et al. 2007) would expect a 
decline in bolometric luminosity at the rate of ^^Co deay, 
0.98 mag (100 day)~^, or faster if the radioactive-decay 
energy is not converted to optical emission with 100% ef- 
ficiency. The shell-shock model (Smith & McCray 2007) 
predicts a rapid decline after ~200 day, though the au- 
thors note that this decline may be offset by the produc- 
tion of a large amount of ^^Ni or continued CSM inter- 
action. The pulsational pair-instability model (Woosley, 
Blinnikov, & Hcgcr 2007) predicts another SN explosion 
at the location of SN 2006gy about 9 years after the ini- 
tial outburst from SN 2006gy. 

More than a year after the explosion. Smith et al. 
(2008b) detected SN 2006gy in the near-infrared (NIR) 
at a luminosity comparable to that of the peak luminos- 
ity of most SNe II. This property had not been predicted 
by any of the models. When coupled with the lack of 
a detection in the radio and X-rays, this led Smith et 
al. (2008b) to conclude that the luminosity could not be 
powered by the continued interaction of the SN ejecta 
with CSM. Kawabata et al. (2009) draw similar conclu- 
sions based on their detection of weak Ha emission at 
a comparable epoch. Smith et al. (2008b) conclude that 
only two possibilities are able to explain the late-time ob- 
servations of SN 2006gy: (i) the explosion produced >2.5 
M0of ^^Ni, possibly as a pair-production SN, which was 
heating dust and consequently generating the large NIR 
excess, or (ii) a massive (~5-10 M©), dusty shell, located 
'--^1 light year from the site of the SN, was being heated 
by the radiation produced at peak, and re-radiated that 
energy as a NIR echo (Dwck 1983). For the first case, if 
the luminosity were powered by radioactive decay, then 
future observations should indicate a continued decline at 
the rate of ^^Co decay. A dust echo, on the other hand, 
would result in a NIR light curve that stays roughly con- 
stant for ~1 2 years before exhibiting a rapid decline. 

Evolved massive stars, such as red super giants and 
luminous blue variables (LBVs), are often observed with 
massive dust shells, so if these stars explode as Type Iln 
SNe one might expect a late-time IR echo. Many SNe Iln 
have been observed to exhibit a late-time NIR excess 
(Gerardy et al. 2002) which in some cases lasted > 1 year. 
This excess has been attributed to NIR echoes (Gerardy 
et al. 2002), though we note that the formation of new 
dust has been argued specifically for SN 1998S (Pozzo 
ct al. 2004) and SN 1995N (Fox et al. 2009). Dusty re- 
gions near the SN should also lead to UV/optical scat- 
tered light echoes (Chevalier 1986); thus, dust is capable 
of providing significant optical and NIR emission at late 
times. The optical decline of Type Iln SNe at late times 
is very heterogeneous (Li et al. 2002), and therefore cau- 



tion must be applied when determining the source of any 
late-time emission. 

In this paper we present new NIR observations of 
SN 2006gy, taken around the peak of optical emission, 
as well as optical and NIR observations obtained more 
than two years after SN 2006gy exploded. In § 2 we 
describe the observations and data reduction. We dis- 
cuss the results in § 3, and in § 4 we offer some con- 
clusions. Throughout this paper we assume that the dis- 
tance to NGC 1260 (the host galaxy of SN 2006gy) is 73. 1 
Mpc, and that the total color excess toward SN 2006gy 
is E{B -V) = 0.72 mag (Smith et al. 2007), of which 
E{B -V) = 0.54 mag is in NGC 1260 and Galactic ex- 
tinction accounts for E{B — V) = 0.18 mag. All spectral 
energy distributions (SEDs) have been corrected for this 
color excess assuming Ry = Av/E{B — V) =3.1 using 
the reddening law of Cardelli, Clayton, & Mathis (1989). 

2. OBSERVATIONS 

NIR observations of SN 2006gy were obtained simulta- 
neously in J, H, and Kg with the Peters Automated In- 
frared Imaging Telescope (PAIRITEL; Bloom ct al. 2006) 
starting on 2006 October 13^ (54 day after explosion^). 
PAIRITEL is a 1.3-m robotic telescope, located on Mt. 
Hopkins, AZ, which obtained observations for 30 min to 
1 hr at each epoch through normal queue-scheduled op- 
erations for the next 24 months. In total there were 124 
epochs obtained. All images were processed via an auto- 
mated pipeline (Bloom et al. 2006). 

Analysis of SN 2006gy has proved challenging because 
the SN is located very close to the nucleus of NGC 1260 
(separation < 1"; Ofek et al. 2007; Smith et al. 2007). 
PAIRITEL has a large native scale of 2" pixel^^, pre- 
cluding spatial resolution of the SN from the nucleus. 
This necessitates image differencing to obtain the light 
curve of the SN. Image subtraction was performed with 
HOTPANTS'^, and the flux in each difference pair was 
determined via aperture photometry at the location of 
the SN. An example subtraction, which clearly shows 
flux from SN 2006gy after the reference image has been 
subtracted, is shown in Figure 1. 

Despite observations extending more than two years 
past the date of discovery, SN 2006gy has not faded be- 
yond the point of detcctability with PAIRITEL. Con- 
sequently, all J, H, and Kg images of SN 2006gy con- 
tain some flux from the SN. We thus adopted the "NN2 
method" of Barris et al. (2005) to determine the relative 
NIR flux changes of the SN. The NN2 method treats 
all images equally and does not require a template im- 
age with no light from the source of interest, SN 2006gy. 
It uses the subtraction of all N{N — 1) possible pairs 
of images to mitigate against possible errors associated 
with the use of a single reference template image. The 
downside to the NN2 method is that it only produces the 
differential flux between each of the N epochs of observa- 
tions. To convert these flux differences into magnitudes 
requires an absolute calibration, which must be obtained 
independently of the results from the NN2 method. Un- 
certainties in the individual subtractions were estimated 

® All dates in this paper are UT unless otherwise noted. 
S Following Smith et al. (2007), we adopt 2006 August 20 as the 
explosion date for SN 2006gy. 

http;//www.astro. washington.edu/users/becker/hotpants. html 
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by measuring the seatter in fake SNe inserted at locations 
having a surface brightness similar to that at SN 2006gy. 

Early-time NIR observations of SN 2006gy (defined 
here as those made before SN 2006gy passed behind the 
Sun during 2007) show a remarkably flat light curve, as 
seen in Figure 2. To transform the relative flux differ- 
ences from the NN2 method to an absolute scale, we sub- 
tracted the archival Two Micron All Sky Survey (2MASS; 
Skrutskie et al. 2006) image from each of the images ob- 
tained on 2006 Nov. 13.24, 14.25, 15.38, and 16.37 to de- 
termine the J, H, and Kg magnitudes of the SN on these 
dates. The mean SN flux was then used to transform the 
relative flux from the NN2 method to the absolute flux of 
the 2MASS system (Cohen, Wheaton, & Megeath 2003). 
PAIRITEL uses the old 2MASS camera and telescope; 
hence, we do not expect any large systematic effects in 
the 2MASS subtractions. The flnal calibrated J, H, and 
Ks photometry is summarized in Table 1. 

On 2006 Nov. 01, Ofek et al. (2007) obtained Palomar 
Hale (5 m) AO images in the J and Ks bands, clearly 
resolving SN 2006gy from the host-galaxy nucleus, and 
mcasiired its flux. In the Ks band our calibration and 
the Ofek et al. (2007) measurement agree to within la, 
while the agreement in the J band is somewhat worse 
(~2c7). Ofek et al. (2007) had only a single 2MASS star 
within the field of view of their AO images, whereas >100 
2MASS stars were used to cahbrate the PAIRITEL im- 
ages; when coupled with the difficulty associated with 
photometry of AO images, this may explain the differ- 
ences between the two measurements. We note that were 
we to adopt the Ofek et al. (2007) measurements as our 
calibration, there would be an overall systematic shift of 
our J and Kg light curves to brighter values, which in 
turn would lead to worse agreement between the NIR 
data and early-time optical spectra (see Figure 4). This 
suggests that our calibration method is sufficient. We 
note that the uncertainties in our photometry are dom- 
inated by the uncertainty in the calibration, which is 
~0.03 mag in J, ~0.06 mag in H, and ~0.04 mag in Kg. 
This uncertainty is the same for all epochs, so a change 
in the calibration would lead to a systematic shift of the 
entire light curve. 

As shown by Smith et al. (2008b), and subsequently 
confirmed by Agnoletto et al. (2009), the NIR evolution 
of SN 2006gy is very slow at late times. Given the rel- 
atively small change in flux, and the reduced signal-to- 
noise ratio following the fading of the SN, we were unable 
to recover reliable flux measurements from PAIRITEL 
data taken after 2007 Sep. Furthermore, unlike the case 
at early times, the SN had faded below the 2MASS detec- 
tion limit, which means that the late-time subtractions 
relative to the 2MASS template image do not yield mean- 
ingful results despite the fact that at Kg « 15 the SN is 
well above the PAIRITEL detection limit. In principle, 
if deep PAIRITEL images are obtained after SN 2006gy 
fades well beyond the detection limit, it should be pos- 
sible to recover the late-time NIR light curve using tem- 
plate images that contain no light from the SN. 

We also observed SN 2006gy with the Near-Infrared 
Camera 2 (NIRC2) using the laser guide star (LGS) AO 
system (Wizinowich et al. 2006) on the 10-m Keck II tele- 
scope in Hawaii on 2008 Aug. 25. We have re-reduced 
the LGS AO images presented by Smith et al. (2008b) 
and derive revised values for the /C'-band magnitudes of 



TABLE 1 

PAIRITEL Observations of SN 2006gy 
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0. 


,09 


13, 


.11 


± 


0, 


.07 


12, 


.93 


± 


0. 


,05 


54093. 


,12 


13. 


,51 


± 


0. 


,11 


13 


.06 


± 


0, 


.06 


12, 


.90 


± 


0. 


,07 


54094. 


.95 


13. 


.68 


± 


0. 


.11 


13, 


.12 


± 


0, 


.07 


12, 


.90 


± 


0. 


.05 


54095. 


.95 


13. 


.66 


± 


0. 


.11 


13. 


.13 


± 


0. 


.07 


13. 


.00 


± 


0. 


.16 


54100. 


.97 


13. 


.79 


± 


0. 


.08 


13 


.15 


± 


0. 


.07 


13. 


.03 


± 


0. 


.13 


54101. 


.42 


13. 


.81 


± 


0. 


.07 


13 


.16 


± 


0. 


.07 


13 


.16 


± 


0. 


.14 



"Midpoint between the first and last exposures in a single stacked 
image. 

''Observed value; not corrected for Galactic extinction. 

SN 2006gy, as summarized in Table 2.^^ We also present 
the flrst measurement of the i?-band flux from the AO 
images taken on 2007 Dec. 2. To obtain this if-band 

measurement, despite a lack of 2MASS stars in the fleld, 
we measured the H — K' color of SN 2006gy relative to 

Latc-timc AO Keck observations were all made in /f'-band. 
The uncertainty associated witli the transformation between K'- 
band and Ks-hanA is small compared to the absolute calibration 
uncertainty and therefore for the late-time AO images K' R! K^. 





Fig. 1. — Example subtraction of PAIRITEL NIR images of SN 2006gy, which have all been registered to the same coordinates. Each 
image is ~2' X 3' in size. Left: PAIRITEL Ks-hand image taken on 2006 Nov. 15. Middle: PAIRITEL Ks-band image taken on 2007 Jan. 
01. Right: Image (left panel) minus reference (middle panel) subtraction image. The image subtraction was performed with HOTPANTS, 
and SN 2006gy is clearly visible in the difference image as a bright point source near the galaxy nucleus. 
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TABLE 3 

HST DAY 825 Observations of SN 2006gy 



50 100 150 200 

Days since explosion [observed frame] 

Fig. 2. — Early photometric evolution of SN 2006gy. Unfiltered 
(approximately R band) observations are taken from Smith et al. 
(2007). PAIRITEL J, H, and Kg observations are from this work. 
We also show the Palomar AO photometry from Ofek et al. (2007). 
The data have not been corrected for extinction in the host or the 
Galaxy. The NIR evolution is remarkably flat, with the R — K 
color increasing steadily after day ~70. This behavior is consistent 
with a cooling blackbody (see text). 



TABLE 2 

Keck AO Observations of SN 2006gy 



date 


epoch^ 


Filter 


mag'' 


(UT) 


(day) 




(Vega) 


2007 Sep. 29 


401 


K' 


14.91 ± 0.17 


2007 Dec. 02 


469 


H 


16.8 ± 0.3 


2007 Dec. 02 


469 


K' 


15.02 ± 0.17 


2008 Aug. 25 


736 


K' 


15.59 ± 0.21 



^Observed-frame day from the adopted explosion date, 2006 Aug. 
20 (Smith et al. 2007). 
'^Observed value; not corrected for Galactic extinction. 

the H — K' color of the host galaxy, and calibrated this 
against the H — Kg color of the galaxy in the archival 
2MASS images. The large uncertainty for this measure- 
ment reflects the accuracy with which we can determine 
the color of the galaxy from the 2MASS images. 

As part of a Hubble Space Telescope (HST) snapshot 
survey (GO-10877; PI Li), SN 2006gy was observed with 
the Wide Field Planetary Camera 2 (WFPC2; Hohzman 
et al. 1995a) using the i^450W^, F555W, F675W, and 
F8UW filters on 2008 Nov. 22. The data were reduced 
in the standard fashion using multidrizzle (Koekemoer 
et al. 2002). NGC 1260 and SN 2006gy were located at 
the center of the PC chip of WFPC2, which has a na- 
tive pixel scale of 0.0455" pixel~^.We follow the recipe of 



Filter 


\ ^ 

■^ccnt 


AA'' 


mag'^ 


Flux'' 




(A) 


(A) 


(Vega) 


(MJy) 


F450W 


4519 


957 


22.16 ± 0.03 


5.97 ± 0.18 


F555W 


5397 


1226 


21.51 ± 0.02 


9.31 ± 0.19 


F675W 


6697 


866 


21.07 ± 0.03 


11.29 ± 0.34 


F8UW 


7924 


1500 


20.91 ± 0.04 


10.80 ± 0.43 



''Central wavelength of the filter, based on WFPC2 calibrations 
presented in Table 8 of Holtzman et al. (1995b). 

''Filter width, based on WFPC2 calibrations presented in Table 8 
of Holtzman et al. (1995b). 

•^Observed value; not corrected for Galactic extinction. 

'^Observed flux using the photometric calibrations presented in 
Table 9 of Holtzman et al. (1995b). 

Dolphin (2000) to do charge transfer efficiency correction 
and photometric reduction of the WFPC2 images. Our 
measurements of the SN 2006gy magnitudes are summa- 
rized in Table 3. A false-color image of this detection is 
shown in Figure 3. 

3. DISCUSSION 

3.1. Early-Time NIR Observations 

While the early evolution of SN 2006gy is remarkably 
flat in the NIR (see Figure 2), we find that these mea- 
surements are consistent with radiation from a cooling 
blackbody. From blackbody flts to optical spectra, Smith 
et al. (2009b) flnd that the temperature of SN 2006gy 
monotonically cools after day ^50, roughly two weeks 
before optical maximum, until day ^^165, where the tem- 
perature levels off at ~6300 K. In the top panel of Fig- 
ure 4 we show the evolution of the SED of SN 2006gy for 
three epochs (day 59, 109, and 135) during our early-time 
NIR observations. For each epoch we show the single- 
component blackbody spectrum, after adopting the tem- 
perature from Smith et al. (2009b) and normalizing the 
spectra to the photometric measurements from Smith et 
al. (2007). Their photometric observations come from a 
series of unflltered observations taken with the Katzman 
Automatic Imaging Telescope (KAIT; Filippenko et al. 
2001), which are best matched by the R band (Riess et 
al. 1999; Li et al. 2003). However, the scatter in the 
transformation between KAIT unfiltered and R can be 
quite large (Ganeshalingam 2009), so we adopt a 0.2 mag 
uncertainty for the calibration of the blackbody spectra. 
We also show the NIR flux during these epochs, deter- 
mined from the absolute calibration of the 2MASS sys- 
tem (Cohen, Wheaton, & Megeath 2003). This shows 
excellent agreement with an extrapolation of the early- 
time spectra of SN 2006gy (Smith et al. 2009b). Between 
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Fig. 3. — Late-time optical and NIR images of SN 2006gy. Top: 
Keck AO H and K' false-color image from Smith et al. {2008b). 
The SN is clearly red in the NIR. Note that the field of view for this 
image is slightly smaller than that of the others. Middle: Keck AO 
K' image of SN 2006gy taken on day 736. The SN is still clearly 
visible in the NIR, and has shown little change since observations 
taken around day 400. Bottom: HST WFPC2 false-color image in- 
cluding the four filters in which we detect SN 2006gy: FibOW (cor- 
responding to blue), i^555iy, F&JbW , and FS14W (corresponding 
to red). The SN is clearly blue in the optical compared to the light 
from the surrounding stars. 

day ~55-135 roughly 2-4% of the bolometric himinosity 
of SN 2006gy was emitted in the NIR. 

In the lower panel of Figure 4 we show the fractional 
excess of the photometric observations relative to the 
single-component blackbodies (note that by definition 
the i?-band excess is set to zero). For clarity we do 
not show the uncertainties associated with the excess, 
but after accounting for the large uncertainty in the R- 
band calibration, each point is within ^la of zero excess. 
Nevertheless, there is an apparent trend that the excess 




4000 6000 8000 10000 20000 



Fig. 4. — SED evolution of SN 2006gy at early times show- 
ing a possible trend towards a growing IR excess above a cooling 
blackbody. Top: the solid, dotted, and dashed lines show a single- 
component blackbody spectrum, normalized to unfiltered observa- 
tions from Smith et al. (2007), on day 59, 109, and 135, respectively. 
Bottom: the fractional excess emission in J, H, and Ks , relative to 
the blackbody model. For clarity we do not include the uncertain- 
ties in these measurements, which are large and dominated by the 
uncertainty in the transformation between unfiltered data and the 
R band (see text). Each point is within r^la of showing no excess; 
however, the NIR excess does appear to grow with time. 

is growing in the H and Kg bands as a function of time. 
The trend may be indicative of radiation from warm dust 
(see below), though we note that this effect would be 
small. If this trend toward a NIR excess is created by 
the same source as the late-time NIR excess, we would 
expect a if^-band excess of ~0.17 mag, comparable to 
the uncertainty in our blackbody calibration. 

3.2. Late-Time NIR Observations 

With data obtained more than a year after explosion, 
Smith et al. (2007) discovered a significant NIR excess 
from SN 2006gy. Our NIR observations taken 736 day af- 
ter explosion show that the A''-band flux from SN 2006gy 
has only faded by a factor of '^2 over the course of the 
previous year. We cannot determine the total IR lu- 
minosity at late times, because our NIR data do not 
cover the peak of the SED. Furthermore, as first noted 
by Smith et al. (2008b), the very red H — K' color at late 
times indicates that the IR emission likely peaks in the 
mid-IR. We can, however, place lower limits on the NIR 
luminosity by assuming that the SED peaks in the K' 
band. Following this assumption, the NIR excess consti- 
tutes a slowly varying luminosity of > 2x10^ Lq for > 

1 yr- 

The most likely explanation for this large luminosity 
in the NIR is warm dust. The observed H — K' color 
from day 436 corresponds to a temperature of ~1000 K, 
assuming a single-temperature blackbody. If the dust is 
radiating as a single blackbody with Tpcak = 1000 K, this 
would increase the above luminosity to > 3 x 10® Lq. 
We note that our NIR observations are only sensitive to 
the warmest dust; it is possible that cooler dust with 
T Ri 600 could dominate the IR emission, in which 
case the luminosity could be significantly higher than the 
values quoted above. 

^ Smith et al. (2008b) show that the combination of the peak 
luminosity and distance to the dust suggest an equilibrium tem- 
perature around 600 K. 
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The location of this dust, and whether it is newly 
formed or pre-existing at the time of the SN explosion, 
remain to be determined. The dust-cooling time is short, 
meaning that a prolonged heat source is needed to ex- 
plain the extended excess. We consider four possibilities 
for heating the dust: (i) radioactive heating from ^^Co 
decay, (ii) coUisional excitation of pre-existing dust, (iii) 
heating via radiation from circumstellar interaction, and 
(iv) a late-time IR dust echo, where pre-existing dust is 
heated by the radiation produced while the SN was near 
its optical peak. 

3.2.1. Radioactive Heating from ^^Co 

Smith et al. (2008b) noted that the observed K'-hmiA 
decay between 2007 September and 2007 December could 
be explained with radioactive heating from a minimum 
of 2.5 Mq of ^^Ni, if a sufficient amount of dust formed in 
order to move the luminosity into the NIR (though they 
favored another interpretation; see below). We show the 
bolometric evolution of SN 2006gy through the first ^800 
day post-explosion, including both optical and NIR de- 
tections, as well as the expected light curve from 2.5 M© 
of ^^Ni, in Figure 5. The early-time measurements (< 
250 day) come from Smith et al. (2009b). The optical 
measurement near day 400 comes from photometric mea- 
surements by Kawabata et al. (2009), while the optical 
luminosity on day 825 comes from a direct integration 
of our HST observations (see § 3.3.1). The three late- 
time NIR luminosities represent lower limits based on 
the measured iiT'-band flux (see above, § 3.2). The decay 
rate of ^^Co, 0.98 mag (100 day)""'^, is much faster than 
the observed decay of the K' flux from SN 2006gy, ~0.2 
mag (100 day)~^. The late-time NIR excess declines at 
a rate that is too slow to be explained by ^^Co heating 
alone. 

The pair- instability model of Nomoto et al. (2007), 
which provided good agreement with the early light curve 
of SN 2006gy after an artificial reduction of the total 
ejecta mass from the evolutionary calculation, was able 
to reproduce the late-time NIR luminosity observed by 
Smith et al. (2008b). This model required less than 
100% efficiency in the conversion of gamma-rays (from 
radioactive decay) to optical/NIR emission, which means 
that the light curve should decay faster than 0.98 mag 
(100 day)-^ The possibility of a pair-instability SN 
was first invoked to explain the large peak luminosity 
of SN 2006gy (Ofek et al. 2007; Smith et al. 2007). This 
scenario would have required the production of > 10 M© 
^^Ni, which in turn would produce a large late-time lu- 
minosity that decays at the rate of ^^Co. The late-time 
NIR luminosity is not accounted for in either the general 
pair-instability models or the artificial model of Nomoto 
et al. (2007); it therefore provides a serious challenge to 
the pair-instability SN hypothesis. ^'^ 

3.2.2. CoUisional Heating of the Dust 

Another possibility is that the dust existed prior to 
the SN explosion, at large distances from the explosion 
site, and was heated via collisions with the expanding 

The pulsational pair-instability model of Woosley, Blinnikov, 
& Heger (2007) has not, however, been excluded; see also Smith et 
al. (2009b). 
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Fig. 5. — Evolution of the bolometric luminosity of SN 2006gy 
during the first 800 day after explosion. Optical data are shown as 
black circles, and NIR data are red squares. Early-time data (< 
250 day) are from Smith et al. (2009b). The optical detection near 
day 400 is from Kawabata et al. (2009). All other data are from 
this work. NIR measurements are only lower limits to the total IR 
luminosity, because our observations do not sample redward of 2.2 
fiUi (see text). The dashed line shows the expected decline if the 
luminosity were powered by 2.5 Mq of ^^Ni. The flat nature of the 
light curve indicates that radioactivity is not the primary source 
of the late-time emission. 

material in the expanding blast wave. The intense ultra- 
violet/optical output from a SN at its peak vaporizes any 
dust in the vicinity of the SN (Dwek 1983). This radia- 
tion near peak creates a dust-free cavity into which the 
SN ejecta may expand at early times; however, the ejecta 
blast wave will eventually reach the edge of the dust-free 
cavity, at which point coUisional excitations may gener- 
ate NIR emission. 

The large peak luminosity of SN 2006gy, 8 x 10^° L© 
(Smith et al. 2009b), provides significant challenges to 
this coUisional heating scenario. Dwek (1983) shows that 
the radius of the dust-free cavity can be determined from 
the peak luminosity and the dust vaporization temper- 
ature, assuming the grain cmissivity is proportional 
to (A/Ao)-i: 



i?i(pc) = 23 



QuLa{rmLQ) 



0.5 



(1) 



where Ri is the radius of the dust-free cavity, Qi, is the 
mean grain cmissivity, Lq is the peak luminosity, and Ty 
is the dust vaporization temperature. Following Dwek 
(1983), if we assume Qi, = 1 and that Aq = 0.2 /im, we 
find that Ri « 1.4 x 10^^ cm, for a vaporization tem- 
perature « 1000 K. If the absorption coefficient is in- 
stead proportional to A~^, then the cavity becomes even 
larger. Based on the observed width of the Ha line from 
SN 2006gy, Smith et al. (2007) estimate the speed of the 
blast wave to be 4000 km s~^, which would mean that 
this material would take ~70 yr to reach the edge of 
the dust-free cavity. Even if there are ejecta traveling 
at more typical SN velocities of 10,000 km s~^, it would 
take over 28 yr for this material to reach the edge of the 
dust-free cavity. Given that the NIR excess is present '^1 
yr following explosion, and that the SN ejecta could not 
reach the dust in such a short time, coUisional heating 
of dust by SN ejecta cannot explain the observed NIR 
excess. 
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3.2.3. Newly Formed Dust 

The emission features that arise in SNe Iln, as in 

SN 2006gy, originate from the interaction of the SN 
ejecta with a dense CSM. Shocked gas, located between 
the forward shock which is plowing into the CSM and 
a reverse shock of the SN ejecta, can cool, and if the 
density is large enough, it may form new dust grains. 
The cooling time for such dust is short (~10-100 s), but 
a persistent heat source can generate an extended pe- 
riod of NIR emission. The physical conditions for post- 
shock dust formation are not necessarily easy to produce: 
strongly interacting SNe are often X-ray sources, which 
could prevent the formation of new dust grains in the 
post-shock gas. In a few cases, however, evidence for the 
formation of dust in the post-shock region has been found 
for SNe 1998S (Pozzo et al. 2004), 2006jc (Smith, Foley, 
& Filippenko 2008), and 2005ip (Smith et al. 2009a; Fox 
et al. 2009), and possibly also SN 2006tf (Smith et al. 
2008a). This newly formed dust is then heated by the 
energy from the shock as it continues to propagate into 
the CSM, which gives rise to the prolonged NIR excess. 

This scenario is difficult to reconcile with the case of 
SN 2006gy, however. Spectra taken around day 200 show 
a decline in the luminosity of the Ha emission line, which 
indicates a reduction in the CSM interaction at this time 
(Smith et al. 2009b). Furthermore, as detailed by Smith 
et al. (2008b), the lack of X-ray, Ha, and radio emission 
a little more than a year after the SN explosion implies 
that the observed NIR luminosity cannot be explained 
by shock- heated dust. In fact, the radio nondctection 
of SN 2006gy continues through day 648 (Bietenholz & 
Bartel 2008). Kawabata et al. (2009) do claim the de- 
tection of very weak Ha emission from SN 2006gy on 
day 394, but the inferred luminosity of ~10^^ erg s~^ 
is nearly three orders of magnitude lower than the ob- 
served NIR luminosity at this time. Kawabata et al. 
(2009) also claim an i?-band detection of 19.4 mag on 
day 394 (though we note that Smith et al. 2008b re- 
port an upper limit of i? > 20 mag, and Agnoletto et 
al. 2009 give an upper limit of i? > 20.3 mag, at a simi- 
lar epoch), which would constitute an optical luminosity 
similar to that observed in the NIR. The source of this 
optical luminosity could potentially be the heat source 
for newly formed dust at day ~400, though it cannot ex- 
plain the similar NIR luminosity seen on day 736. While 
SN 2006gy was detected in the optical at > 800 day to 
have a luminosity similar to that seen in the NIR, the 
very blue nature of optical data suggest that it is a scat- 
tered light echo (see below) of UV/optical emission from 
the SN peak. The light echo must originate outside the 
dust-forming region and therefore cannot heat the dust. 
While we cannot strongly rule out dust formation in the 
post-shock region of SN 2006gy it is implausible to ex- 
plain the late- time NIR luminosity via new dust, because 
there are no signs of a heating source that lasts for > 1 

yr- 

3.2.4. NIR Dust Echo 

Perhaps the most natural explanation for the late-time 
NIR excess is a dust-heated echo, as first discussed by 
Smith et al. (2008b). The qualitative behavior of IR 
echoes predicted by Dwek (1983) provides an excellent 
match to the observed behavior of SN 2006gy. The IR 



emission comes from dust near the explosion site which 
is heated by the early-time UV/optical emission from the 
SN. The model predicts a fast rise in the IR, followed by 
an extended plateau (see below) before the IR luminosity 
follows a rapid decline. SN 2006gy exhibits a very slow 
decline in the NIR, only a factor of ~2 in flux over the 
coiirse of ^1 yr, consistent with being on the extended 
plateau phase as predicted by the model. 

The plateau occurs because the hottest dust, which 
exists right at the edge of the dust-free cavity described 
above, dominates the emission. According to the view 
of a distant observer, the emitting volume is a series of 
paraboloid light fronts that expand throughout the dust- 
free cavity (see Figure 1 of Dwek 1983 or Figure 4 of 
Smith et al. 2008b). As a paraboloid expands it contin- 
ually heats dust at Ri, which is what gives rise to the 
plateau in the light curve. After a time ^2Ri/c, the ver- 
tex of a given paraboloid will reach the back edge of the 
cavity and will no longer be heating dust at Ri. Emis- 
sion from the dust will be dominated by the UV/optical 
radiation produced by the SN at peak, and once this ra- 
diation sweeps past the back edge of the cavity the IR 
luminosity begins to rapidly decline. Therefore, if we 
know the duration of the plateau phase we can deter- 
mine the size of the dust-free cavity. Constraints on the 
cavity size can be determined by assuming that the echo 
starts at our first K' observation on day 405 and ends 
on our last observation on day 736. In this case we find 
that i?i w 4 X 10^'' cm, however, it would be contrived if 
our observations were to perfectly bookend the plateau 
phase of the echo: the echo almost certainly started prior 
to day 405 and continued beyond day 736. If we instead 
assume that the slight NIR excess seen at day ~100 (see 
Figure 2) is the rise of an IR echo, then we find that 
iii i=s 8 X 10^^ cm, which shows reasonable agreement 
with the value we calculated for the dust-free cavity in 
§ 3.2.2, Ri « 1.4 X 10"'^* cm, considering the uncertain- 
ties in the dust temperature and the likelihood that the 
plateau extends beyond day 736. 

Our NIR observations are only sensitive to the warmest 
dust, but if we assume that the NIR luminosity peaks 
in the K' band, then a lower limit to the total energy 
emitted in the IR is E^i « 2 x 10^ Lq x 600 d, or > 
4x10^^ erg. This is comparable to the canonical optical 
output of a normal SN II. Typically, this would pose 
a significant problem for the IR echo hypothesis as it 
would require order unity of the original optical radiation 
from the SN to be absorbed and re-radiated in the NIR. 
Indeed, this is precisely the argument used by Fox ct al. 
(2009) to argue against the light-echo hypothesis for SN 
2005ip. In the case of SN 2006gy this is not a problem, 
however, because the NIR emission falls well within the 
available energy budget to heat the dust. After applying 
a bolometric correction, Smith et al. (2009b) find that 

2.5 X 10^^ erg were emitted by SN 2006gy during the 
first ~220 day of the light-curve evolution, which means 
that less than 2% of this energy has been absorbed and 
re-radiated by the dust. This value is comparable to 
those found by Dwek (1983) for SNe 1979C and 1980K. 

If we assume that the dust radiates as a blackbody, 
we can estimate the optical depth of the dust shell as 
T EiR/{E+Em) « 0.016. We note that like Em above, 
this value constitutes a lower limit for r. Even if the NIR 
luminosity is a factor of 10 larger than our lower limits, 
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this would imply r w 0.1, which corresponds to an optical 
extinction of Ay — l.OSGr « 0.1. This is far less than 
the total observed extinction from the host, Ay^host ~ 
1.67 (Smith et al. 2007), suggesting that the dusty shell 
responsible for the IR echo cannot produce all of the non- 
Galactic reddening observed toward SN 2006gy. That 
there is a considerable amount of dust along the line of 
sight, which is not in the GSM of SN 2006gy, should not 
come as a surprise given that there is a prominent dust 
lane in NGG 1260 (see Figure 3). 

With this information in hand, it is now possible to 
estimate the total dust mass of the shell that creates the 
NIR echo; according to Dwek (1983), it is simply related 
to the dust properties, the optical depth of the dust, and 
the inner and outer radii of the shell. Above we estimate 
the inner radius of the dust shell to be '-^8x10^^ cm. 
The value of the outer radius can be deciphered from 
the decay of the NIR light curve following the plateau. 
As we have not yet observed the NIR decay, we adopt 
an outer radius of 2Ri, which minimizes the dust mass 
in the shell. Finally, we follow the same assumptions as 
Dwek (1983) about the dust properties: a grain density 
Pgr = 3 g cm~^, grain radius a = 0.1/im, and mean ab- 
sorption efficiency Q^, = 1. With these asumptions we 
find a total dust mass of Md ~ 0.1 Mq. Assuming a 
typical value for the dust-to-gas ratio, 1:100, the total 
mass of the circumstellar shell is ~10 Mq. This value 
shows good agreement with the initial estimate of 5-10 
M© from Smith et al. (2008b), which was based on only 
observations taken between 400 and 465 day after explo- 
sion. We note that ~10 Mq constitutes a lower limit to 
the total mass of the shell, because the actual value of 
the duration of the plateau and E'ir could potentially be 
much larger than the values we adopted above. 

At a distance of ^lO""^* cm from the SN explosion site, 
it is not immediately evident that the dust giving rise to 
the NIR echo is part of the progenitor's GSM or the local 
interstellar medium (ISM). Smith et al. (2008b) argue 
that the large dust mass needed at ~10^* cm could be 
explained if the progenitor passed through a phase of 
eruptive mass loss ^^1000-1500 yr prior to explosion, in 
an event similar to observed outbursts from rj Car. Pre- 
existing dust in the walls of a giant HII region, where 
a massive star like SN2006gy's progenitor might have 
lived, cannot account for the IR echo because at typical 
distances of ~10 pc or more from the SN, the dust would 
be far too cold to reproduce the NIR excess. 

3.3. Observed Color Evolution from HST Observations 

Ghevalier (1986) predicts that any SN with an IR echo 
due to pre-existing GSM dust should also show a faint 
scattered light echo in the optical as well. He shows that 
this dust can lead to optical emission characterized by a 
A^" scattering law, where a is some value between 1 and 
2. More than a year after explosion SN 2006gy had faded 
rapidly in the optical (Smith et al. 2008b; Agnoletto et 
al. 2009; Kawabata et al. 2009). Therefore, we obtained 
HST images to search for a possible late-time scattered 
light echo. 

Our optical detection on day 825 shows a significant 
change in the optical decline rate of SN 2006gy. Between 
day 200-400 the SN faded by > 3 mag in the optical 
(Smith et al. 2008b), whereas the decline between day 
400-825 was only ^^1 mag. In addition to this reduced 
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Fig. 6. — Late-time SED of SN 2006gyshowing evidence for 
two distinct emission components. We show the flux in the optical 
{HST, day 825) and NIR (Keck AO, day 469 and 736), corrected for 
host-galaxy and Milky Way reddening. The solid blue line shows a 
12,000 K blackbody, normalized to the F8UW detection. The red 
dotted line shows a single-component blackbody fit to the H and 
K'-band detections from day 469 (note that this fit has zero degrees 
of freedom). The IR evolution is slow, and the detection from day 
736 indicates the degree of fading in the NIR. The SED shows a 
rapid decline toward the red portion of the optical, as well as a 
strong rise toward the red within the NIR. These characteristics 
are precisely those predicted by the IR and scattered light echo 
models. 

rate of decline, the SN underwent significant color evolu- 
tion. Assuming that F555W^ and F814W^ approximate 
the V and / bands, respectively, we find that V—I — 0.60 
mag on day 825. By contrast, with a 12000 K color tem- 
perature at peak (Smith et al. 2009b), the observe color 
of SN 2006gy would have been F - / w 1.01 mag. Thus, 
the spectrum on day 825 is bluer than the spectrum at 
peak. Indeed, the '^0.4 mag color change shows excel- 
lent agreement with a scattering law. A scatter- 
ing law, on the other hand, would predict a —0.83 mag 
change to the V — I color, or roughly V — I ^ 0.20 mag 
during the late-time epochs. 

This behavior is similar to what one would expect from 
a scattered light echo: as in the case of an IR echo, when 
successive paraboloids sweep out to progressively larger 
radii, dust in the circumstellar (or in some cases the in- 
terstellar) environment can scatter that light toward the 
observer, creating a plateau in the optical light curve of 
the SN. At the same time, scattering preferentially se- 
lects shorter wavelengths, which results in a spectrum 
that is bluer than that of the SN at peak. This is the 
precise behavior observed for the Type la SNe 1991T 
(Schmidt et al. 1994), 1998bu (Gappellaro et al. 2001), 
and 2006X (Wang et al. 2008) which both showed signif- 
icant departures from the expected decline rate of Type 
a SNe. Spectra taken roughly two years after maximum, 
for SNe 1991T and 1998bu, and -10 months after max- 
imum for SN 2006X, show emission features similar to 
those seen in the spectra near peak superposed on top of 
a blue continuum for both these SNe. Light echoes have 
also been observed around a number of core-collapse SNe, 
e.g. SN 2003gd (Sugerman 2005; Van Dyk, Li, & Fihp- 
penko 2006) and SN 1993 J (Sugerman & Grotts 2002), 
though we note that late-time spectra of the light echoes 
do not exist in these cases. 

In Figure 6 we show the late-time SED of SN 2006gy, 
including our NIR detections on day 469 and 736 and 
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Fig. 7. — Evidence for a scattered light echo from the blue optical 
emission from the SN. The _f/5T/WFPC2 photometry of SN 2006gy 
(black dots) are compared to observed and reficeted spectra, both 
normahzed to the HST F&IAW flux density. The gray line is the 
observed spectrum of SN 2006gy on day 71 during its peak lumi- 
nosity phase, from Smith et al. (2009b). We plot the observed 
wavelength (i.e., not corrected for the SN redshift), and the spec- 
trum has no correction for reddening. The black line (and dotted 
lines above and below) show the same spectrum, but adopting an 
assumed wavelength dependence for the dust reflectcUice propor- 
tional to A-l-2(±015). 

the optical detections from day 825. The solid blue line 
represents a single-component blackbody at T = 12,000 
K, roughly the peak temperature of SN 2006gy (Smith 
ct al. 2009b), normalized to our F814VK detection. The 
SED clearly shows that the late-time emission is signifi- 
cantly bluer than a blackbody spectrum. We note that 
the early-time spectra of SN 2006gy show considerable 
line blanketing blueward of ~4000 A, such that the emis- 
sion decreases blueward of our _F450T4^ detection. The 
NIR emission rapidly rises toward the IR, suggesting that 
there is considerable IR emission to which our NIR mea- 
surements are not sensitive. 

Figure 7 illustrates the observed HST photometry 
along with the spectrum of SN 2006gy obtained near 
peak (Smith et al. 2009b). Note that no correction 
for reddening has been made. The photometry shows 
excellent agreement with a A~" scattering law, where 
a = 1.2 ±0.15. The lack of an actual reflection spec- 
trum precludes a more detailed model of the scattering 
dust at this time. Given that the other typical scenarios 
for a large late-time luminosity from a SN, ^^Co decay 
and extended CSM interaction, have been ruled out, we 
conclude that the late-time HST detection is a scattered 
light echo. 

3.4. Location of the Scattering Dust 

The location of the scattering dust is of interest as it 
can help determine properties of the environment local to 
SN 2006gy. Above we show that the CSM dust responsi- 
ble for the NIR echo cannot account for all the reddening 
observed toward SN 2006gy, suggesting that the SN light 
passes through external, non-CSM dusty regions in the 
host galaxy prior to reaching us. Therefore, the scat- 
tered light echo could originate in a number of different 



locations. 

With only a single epoch of the unresolved scattered 
light echo, a precise determination of the location of the 
scattering dust cannot be obtained. Nevertheless, we can 
gain some insight into the dust location based on the 
total optical and NIR emission. For the case where the 
NIR-echo dust also gives rise to the scattered light echo, 
and assuming the SN can be modeled as a single short 
pulse of UV/optical radiation, the ratio of the scattered 
light luminosity (Lg) to the IR luminosity at late- times 
can be shown to iDe (Chevalier 1986): 

Ls _ LP A {ctf\ {1-9? 

/ (l-g)2 + 2g {l-gf+gct/R2 \ 

1 [(1 - gf + 45] 1/2 [(1 - g)2 + 2gct/R2Y/^ J ' 

where uj is the dust albedo, c is the speed of light, t is the 
time since the SN explosion, R2 is outer radius of the dust 
shell, and 5 is a measure of the degree of forward scatter- 
ing. The case 5 = corresponds to isotropic scattering 
and g = 1 is completely forward scattering. Note that 
in Equation 2 wc have assumed that ct/2 < Ri, where 
Ri is the inner radius of the dust shell, and that the 
albedo is independent of frequency. Empirical estimates 
and numerical calculations (White 1979) show g « 0.6. 
For the dust albedo we adopt ui « 0.6 (Mathis, Rumpl, 
& Nordsieck 1977). As mentioned above, we cannot con- 
strain i?2 with our current observations, but we do know 
that i?2 must be greater than Ri , and the value of Equa- 
tion 2 does not change when R2 00. For reasonable 
limits to the outer radius, 1.5i?i < R2 < 00, we find that 
Ls/Lin « 0.88 0.60. 

Direct integration of the HST optical flux yields an 
optical luminosity of ~(7.4±0.5)xlO'' Lq. This value is 
somewhat uncertain because our data do not extend into 
the ultraviolet; however, optical spectra taken near peak 
show that line blanketing severely reduces the flux blue- 
ward of ~4000 A (Smith et al. 2009b), so this uncertainty 
should not have a significant effect on the total scattered 
light luminosity. If we assume that the NIR flux contin- 
ues to decline at the same rate as that observed between 
day 405 and 736, then we find a lower limit to the IR lu- 
minosity (again because our NIR observations only probe 
emission from the hottest dust) of Lir > 8.6 x lO'' Lq 
on day 825. Therefore, on day 825, Ls/Lm < 0.86, con- 
sistent with the predicted values above in the case where 
the scattering dust is the same dust responsible for the 
IR echo. 

The alternative possibility is that the scattering dust 
is not located in the CSM, instead existing at some other 
location in NGC 1260 along the line of sight. Cappellaro 
ct al. (2001) inferred that this was the case with SNe 
1991T and 1998bu, and like those two SNe we know that 
SN 2006gy has a significant amount of dust located along 
the line of sight. As shown by Cappellaro et al. (2001), if 
we approximate the SN light curve as a short pulse with 
duration AtsN, then the scattered echo luminosity is 

LccUt) = LsNAfsN,/(t), (3) 
where LsnA^sn can be obtained via direct integration 
of the observed light curve, Lgu^tsu = Jf^°° LsN{t)dt, 
and f{t) is the fraction of light scattered to the observer. 
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Under the assumption that light is being scattered by 
dust in a sheet of thickness A_D <C -D, where D is the 
distance between the dust and the SN, then (Chevalier 
1986; Cappellaro et al. 2001) 



8nD + ct{i+g2 _2g[D/{D + ct)]f/^' 

where r is the optical depth of the dust. Using the rela- 
tion Av = 1.086ry (Cox 2000), and the fact that Ay out- 
side the dusty CSM is ^1.6 mag, we find t w 1.5. From 
Smith et al. (2009b) we know that isNAtsN « 2.5 x 10^^ 
erg. Therefore, using Lecho ~ 2.9 x lO"'^ erg s~^, we find 
(from Equation 4) that dust located at Z? « 20 pc can 
explain the observed scattered light luminosity. 

Our observations appear to be consistent with one of 
two scenarios: either (i) the scattered light echo is due 
to CSM dust, which has also been heated and is radiat- 
ing in the IR, or (ii) a dusty region ~20 pc from the SN 
is responsible for the late-time optical luminosity. This 
latter scenario seems very plausible if the progenitor of 
SN 2006gy was a very massive star, > 100 (Smith et 
al. 2009b). If this were the case we might expect that 
the progenitor resides in the center of a giant HII region 
with multiple dense, dusty clouds nearby. This scenario 
might be very similar to the Carina nebula surrounding rj 
Car, which has multiple dusty molecular clouds only 10- 
20 pc from rj Car (e.g. Smith & Brooks 2007). Perhaps 
given a similar progenitor birth environment, evidence 
for patchy dust at the tens of pc scale around GRB en- 
vironments also appears to be secure (e.g. Perley et al. 
2009). If the progenitor of SN 2006gy were surromided 
by such a nebula then our assumption of a single, thin 
scattering surface would no longer be valid, and the scat- 
tering dust could potentially be located at a number of 
different locations. 

With only a single epoch of observations, we cannot 
distinguish between the two possible dust locations. One 
additional epoch of optical imaging should prove suffi- 
cient to determine which scenario is correct. As long as 
D ct, then f{t) is roughly constant, whereas the mod- 
els of Chevalier (1986) show a continuous decline in the 
scattered light echoes when g = 0.6 and the scattering 
dust is in the CSM. A relatively flat optical light curve 
would put the dust at ~20 pc, whereas significant decline 
from the observed day 825 flux would suggest CSM dust. 
Finally, we note that the prospects for fully resolving the 
SN 2006gy light echo are dim: at a distance of 73 Mpc it 
would take several decades before the light echo could be 
resolved with high spatial resolution images (i.e., compa- 
rable to HST in the optical). 

4. CONCLUSIONS AND FUTURE WORK 

We have presented new observations of SN 2006gy, in- 
cluding early- time NIR data from PAIRITEL, as well as 
optical and NIR detections more than two years after 
the explosion. The early NIR data show good agreement 
with a single-component, cooling blackbody between day 
54 and 133. There is tentative evidence for a slight NIR 
excess toward the end of this period. The late-time NIR 
and optical detections show a rate of decline which is 
slower than that of ^^Co, indicating that the luminosity 
is not generated via heating by radioactive decay. This 
slow rate of decline is not expected in the pair-instability 



SN models. Rather, the NIR observations can best be 
explained with an IR echo, and the optical detections 
are most compatible with a scattered light echo. Both 
IR and scattered light echoes should be present if there 
is a massive dust shell external to the SN (Dwek 1983; 
Chevalier 1986). 

At the time of discovery SN 2006gy was the most 
luminous SN known (Ofek et al. 2007; Smith et al. 
2007), which when combined with its proximity to the 
host-galaxy nucleus (^300 pc), suggested that perhaps 
there was something unique about star formation and 
SNe in that region. Since the discovery of SN 2006gy 
several other VLSNe have been found, including SNe 
2005ap (Quimby et al. 2007), 2006tf (Smith et al. 2008a), 
2008am (Yuan et al. 2008), and 2008es (Miller et al. 2009; 
Gezari et al. 2009). In each of these cases the SN has been 
localized well outside the central kpc of the host galaxy, 
with the exception of SN 2008es, for which a definitive 
host has not yet been identified (Miller et al. 2009; Gezari 
et al. 2009). This suggests that the location of SN 2006gy 
was not in fact special. Current and future surveys will 
dramatically increase the rate of discovery of VLSNe, 
which should, in turn, allow for a more systematic study 
of their host environments (for a discussion of the rates 
of VLSNe see Miller et al. 2009). 

Given our interpretation for the late-time optical and 
NIR emission arising from a dust shell at ^8x10^^ cm 
from the SN, what might we expect from future observa- 
tions of SN 2006gy? The IR-echo models of Dwek (1983) 
show that after a time t « 2Ri/c, the IR luminosity falls 
off the plateau and exhibits significant decline. Based on 
our predicted value for the size of the dust-free cavity (see 
§3.2.2), we would expect that the plateau phase of the IR 
echo has ended, and that the IR flux is now in steady de- 
cline. Future observations of this decline will place limits 
on the outer extent of the dust shell, which will in turn 
provide better limits on the total mass of the dust shell. 
Our ability to predict the behavior of the scattered light 
echo is limited by the current degeneracy concerning the 
physical location of the scattering dust. If the same dust 
is responsible for both the IR and scattered light echoes, 
then we would expect that the optical emission is now 
fading, in a manner similar to what ought to be seen in 
the NIR. Our observations are also consistent with dust 
located at ~20 pc from the SN, in which case we might 
expect the scattered light echo to continue at roughly a 
constant flux for several more years. Continued obser- 
vations with HST would allow us to distinguish between 
these two situations. 

The shell-shock model of Smith & McCray (2007) may 
also explain the peak luminosity of several other VLSNe: 
2005ap, 2006tf, 2008es (Quimby et al. 2007; Smith et 
al. 2008a; Miller et al. 2009; Smith et al. 2009b). This 
model requires large CSM densities at a distance of a few 
X 10^^ cm, which can be accomplished if the SN progen- 
itor undergoes significant (~0.1-1 Mq yr~^) mass loss 
during the few decades prior to explosion (Smith et al. 
2007). The most luminous LBVs have observed mass 
sheUs with >10 Mq (Smith & Owocki 2006 and refer- 
ences therein) , indicative of giant eruptions (Humphreys 
& Davidson 1994) but not LBVs in their more typical 
state with less violent wind variability (see Smith, Vink, 
& de Koter 2004 for a general reference on LBV winds). 
The possible connection between large CSM densities 
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and VLSNe may suggest that their progenitors are LBVs 
(Smith et al. 2009b). For the case of SN 2006gy, Smith et 
aL (2008b) argue that the dusty shell at ~10^* cm could 
exist if the progenitor underwent an LBV-like eruptive 
mass-loss phase ~1200 yr prior to the SN. IR echoes for 
other VLSNe have not been reported to date, but we 
strongly encourage a search for them. If a substantial 
fraction of them exhibit late-time characteristics simi- 
lar to SN 2006gy, this may be suggestive of a common 
timescale, ~1000 years and -^10 years, for extreme mass 
loss in the progenitors of VLSNe. The combination of 
a VLSNe and IR echo would point to multiple phases of 
eruptive mass loss, which is reminiscent of LBV behavior. 
The instability driving these eruptions still remains un- 
clear (one possibility is the pulsation pair-instability de- 
scribed by Woosley, Blinnikov, & Heger 2007 to explain 
SN 2006gy), and we encourage future observational and 
theoretical work to better characterize these systems. 
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